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Statistical properties of effective mass are analyzed. We show from a general ground that effective mass as a 
function of time should not exhibit long plateaux whatever high statistics simulations are made: the mass should 
fluctuate beyond the one standard deviation of error bars after a few time slices for large times where the ground 
state dominates. This explains the difficulty of obtaining long plateaux experienced in previous simulations. 
Implications of the observation for global fits are discussed, and results for hadron masses are presented. 



1. Introduction 

Calculation of hadron masses constitutes a ba- 
sic part in virtually all problems of lattice QCD 
simulations. A quantity always examined in such 
calculations is the effective mass mcff{t). The ex- 
istence of a plateau in mcs{t) as a function of time 
t is regarded as a dual measure for the statistical 
quality of data and the minimum time separa- 
tion beyond which the ground state dominates. 
In this way the plateau is used as a guide in the 
choice of the time interval for a global fitting 
to extract hadron masses from propagators. In 
practice, however, a long plateau has rarely been 
seen. Even in the best previous efforts toward 
high statistics simulations effective masses, 
particularly for p meson and the nucleon, almost 
always deviate from a plateau beyond error bars 
after 5 or 6 times slices. Ensuing uncertainties 
in the choice of the fitting range and fitted val- 
ues of hadron masses have represented a severe 
hindrance factor in attempts toward high preci- 
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sion determination of hadron masses, especially 
for light hadronsQ. 

In this report we present an analysis on the 
origin of this problem. Our study is based on 
hadron propagators for the Wilson quark action 
at K = 0.1545,0.1550,0.1555 evaluated on 1000 
quenched gauge configurations on a 24^ X 64 lat- 
tice at /? = 6.0, which were generated with the 
5-hit pseudoheatbath algorithm at 2000 sweep in- 
tervals. The central value K — 0.1550 has been 
chosen to facilitate a comparison with previous 
high statistics studies which employ up to 

400 configurations for the same spatial sizeQ or 
a larger size of 32'^ . 

This calculation is one of the first QCD runs 
carried out by the JLQCD Collaboration on 
VPP500/80 at KEK which started operation in 
January 1995. The machine consists of 80 pro- 
cessing elements (PE's), each with the peak speed 
of 1.6GFL0PS and 256MBytes of memory, con- 
nected by a crossbar switch. Our run used 64 
processors, on which our code for heatbath and 
red/black minimal residual solver sustained the 
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Figure 1. Effective mass for (a) tt, (b) p and (c) N at P — 6.0 and K = 0.1550 on a 24^^ x 64 lattice with 
1000 configurations. Errors are estimated by a single elimination jackknife procedure. 



speed of 1.0 - 1.3 GFLOPS/PE. The configura- 
tion generation of 2000 sweeps and a calculation 
of standard meson and baryon propagators for the 
point and wall sources on the final configuration 
was made in about 30 minutes, so that the entire 
run took about 20 days to complete. 

2. Analysis of effective mass 

In Fig. 1^ we show the effective mass for tt, p and 
nucleon (N) obtained with the point (open circles) 
or the wall (filled circles) source dX K — 0.1550 
or m-^/mp « 0.7, which roughly corresponds to 
strange quark mass. A striking feature observed 
in these plots is the presence of fluctuations which 
exceed the level of one standard deviation after a 
few time slices even for tt and with 1000 configu- 
rations. 

In fact it is possible to understand this be- 
havior by a simple statistical analysis. To show 
this, let G{t) be the average for a hadron prop- 
agator over N independent configurations, and 
Gtruo(i) = {G{t)) be the true propagator. Ac- 
cording to the central limit theorem, the differ- 
ence 5G{t) = G[t) — Gtruo(^) obeys the distribu- 
tion. 



P[G] cx exp 



N 



-Y,5G{t)C-\t,t')5G{t') 



,(1) 



where the covariance matrix C is defined by 
C{t,t') = {G{t)G{t'))-{G{t)){G{t')). 



(2) 



Let us introduce a normalized covariance matrix 
C{t,t') 



Cit,t') = 



^C{t,t)C{t',t') 



(3) 



and denote the eigenvalues and normalized eigen- 
vectors of C{t, t') by Xi and ei{t) (i = 1, • • • , Lt) 
with Lt the temporal lattice size. An eigenvector 
decomposition of the averaged propagator leads 
to the formula. 



G{t) = Gtrueit) + 



C{t,t ) ^ 

N 



i=l 



(4) 

Lt) are in- 



where, according to (|lj), (i = 1, 
dependent Gaussian random numbers. 

Given a model of the true propagator Gtiuc(0 
and a measured value of the covariance matrix, 
this formula allows us to generate "simulated" 
samples of the averaged propagator G{t) by gen- 
erating a set of Lt Gaussian random numbers . 

For the model of Gtrue(0 take a double hy- 
perbolic cosine form. 



Gtruc(t) - ^e-"* 



Z'e- 



+ {t^Lt-t), (5) 



where the masses and residues are determined 
by a fit of the measured propagator over 
Q < t < Lt/2. In Fig. || we show four exam- 
ples of effective mass for tt and p for the point 
source calculated from "simulated" propagators 
(solid curves). Dotted lines represent values for 
the model (||). These examples clearly demon- 
strate that fluctuations as observed in Fig. |l| are 
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Figure 2. Four examples of effective mass for tt (top row) and p (bottom row) for point source calculated 
from "simulated" propagators (solid curves) as compared to the measured values for K = 0.1550 (shown 
with error bars) taken from Fig. |^. Dotted lines correspond to the model (H). 



a typical occurrence. Indeed some of the exam- 
ples are amusingly similar to the measured value. 

Let us add a remark that the diagonal of 
the covariance matrix is expected to behave as 
C(i, t) cx exp(— at) with a = 2mT^ for n and p and 
a — Sm^r for with which our data are con- 

sistent. Combined with (^) this explains why the 
magnitude of fluctuation of increases rapidly 
for p and N toward large times, while it stays 
roughly constant for tt. 

We now reexamine the concept of a plateau in 
the light of the above analysis. Let us define a 
plateau over the time interval ti < t < t2 hy 
the condition that TOefi(i) over this interval falls 
within a band of width 2(5rn centered at m, i.e., 
m — 5m < moff(t) < m + dm for t = ^i, • • • , 
To calculate the probability for the occurrence 
of such a plateau, we make a change of variable 
P[G] P[m cs{ti) , ■ ■ ■ ,mctf {t2)] in (|l|), defining 
meff(i) = logG{t)/G{t + l) (effects of the peri- 
odic boundary condition are small in the numer- 
ical results below). The probability is then given 
by an integral of P[mt-^ , ■ • ■ , ^12] over m — 5m < 
mti < m + 5m{i = ti, • • • , 12). 

In Fig. H we plot the probability of finding a 
plateau of length tp{— t2 ~ H + \) for tt and p 
calculated with the measured covariance matrix 
for K = 0.1550 for the case of the wall source. 
We fix t\ = 15, since we have chosen the ground 
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Figure 3. Probability for plateau of length tp 
starting at ti — 15 for tt and p. Open and filled 
symbols correspond to the choice of 5m equal to 
one or two standard deviations of rrioff (^i = 15). 

state mass m by a fit with a single hyperbolic 
cosine over lb <t < Lt/2. For 5m we choose one 
(open symbols) or two (filled symbols) standard 
deviation of mcs_{t) at t = ti. Values for p for 
t > 7 are not shown because of the poor quality 
of the covariance matrix for large times. 

This figure shows that one .should not expect 
to see a plateau at the level of one standard de- 
viation even for tt: the probability drops below 
50% after 2 time slices. Allowing for a deviation 
of two standard deviations, a plateau of length 
10 for TT (3 ~ 4 units in terms of the correlation 
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Table 1 

Hadron masses in lattice units at = 6.0. For our 
results and those of QCDPAX the first and second 
row correspond to point and wall source. APE 
used multi-origin 1^ source, while LANL results 
are combined estimates from wall and Wuppertal- 
smeared source. 

K TT ~p N 

0.1545 0.33075(40) 0.4441(6) 0.6828(21) 

0.33076(28) 0.4425(10) 0.6777(21) 

0.1550 0.29642(42) 0.4231(16) 0.6451(11) 

0.29642(27) 0.4220(12) 0.6393(27) 

0.1555 0.25867(46) 0.4019(20) 0.6066(13) 

0.25864(33) 0.4016(17) 0.6003(37) 

Previous results at = 0.1550 

^PE 0.298(2) 0.429(3) 0.647(6) 

QCDPAX 0.2960(8) 0.4201(29) 0.6403(50) 

0.2964(6) 0.4228(19) 0.6307(39) 

LANL 0.297(1) 0.422(3) 0.641(4) 

length = l/mTr) becomes probable, while for p 
the probability decreases to 60 — 50% already at 
t « 6 — 7 (i/^p ~ 3). These features are consistent 
with actual examples of effective mass shown in 
Fig. |l|, and also with those of previous high statis- 
tics simulations 

We should emphasize that the pattern of fluctu- 
ations of effective mass does not change when one 
increases the number of configurations N , except 
that the magnitude scales down as 1/ y/N. At the 
same time statistical errors estimated for meff also 
decreases as 1/ ^/N . In this sense higher statistics 
does not lead to a longer or better plateau. In 
other words we can ask for such an improvement 
only within a fixed magnitude of the absolute er- 
ror {e.g., 1% of mass). Let us add that the use 
of larger spatial volumes and improved operators 
having a larger coupling to the ground state also 
help to obtain a better plateau only in the latter 
sense. 

3. fits for hadron masses 

Our analysis should have made it clear that re- 
stricting the fitting range of a global fit to the 
time interval of an apparent plateau is not well 
founded. In fact if one repeats a simulation with 



a different sequence of random numbers the effec- 
tive mass will generally exhibit a plateau at some 
other time interval at a different value. This indi- 
cates that it is more reasonable to take the mini- 
mum time tmm at the time slice where the dom- 
inance of the ground state is reasonably ensured 
{e.g., by the overlap of mries for the point and wall 
sources), and to extend the fitting range to large 
times as long as statistical fluctuations do not be- 
come unacceptably large, without resorting to the 
presence of plateau. The value of correlated 
should tell whether the choice is reasonable. Of 
course the fitted values of hadron masses vary de- 
pending on the choice of the interval. However, 
this is an uncertainty which can only be reduced 
by an improved measurement of hadron propaga- 
tors. 

In Table we present hadron masses obtained 
by a correlated fit over the interval 15 < 
t < Lt/2 = 32 with a single hyperbolic cosine 
for TT and p and with a single exponential for 
N . Errors correspond to an increase of by 
one. For tt and p our results obtained with the 
point and wall sources are mutually in agree- 
ment. For K = 0.1550 previous results from 
the QCDPAX Collaboration!^ and from the Los 
Alamos group[Q are consistent with ours, while 
those from APE are higher. The case of nucleon 
is problematical. We find that the fit is not 
very stable and that the point and wall results do 
not agree within the error. Effort toward improv- 
ing baryon operators will be needed for a precise 
determination of baryon masses even in the region 
of strange quark. 
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